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Effect of current hysteresis on the spin polarization of current in a paramagnetic 

resonant tunnehng diode 

P. Wojcik, B.J. SpisakQ M. Woloszyn, and J. Adamowski 

Faculty of Physics and Applied Computer Science, 
AGH University of Science and Technology, al. Mickiewicza 30, 30-059 Krakow, Poland 

A spin-dependent quantum transport is investigated in a paramagnetic resonant tunneling diode 
(RTD) based on a Zni_a;Mna;Se/ZnBeSe heterostructure. Using the Wigner-Poisson method and 
assuming the two-current model we have calculated the current-voltage characteristics, potential 
energy profiles and electron density distributions for spin-up and spin-down electron current in an 
external magnetic field. We have found that - for both the spin-polarized currents - two types of the 
current hysteresis appear on the current-voltage characteristics. The current hysteresis of the first 
type occurs at the bias voltage below the resonant current peak and results from the accumulation 
of electrons in the quantum well layer. The current hysteresis of the second type appears at the 
bias voltage above the resonant current peak and is caused by the creation of the quasi-bound state 
in the left contact region and the resonant tunneling through this quasi-bound state. The physical 
interpretation of both the types of the current hysteresis is further supported by the analysis of the 
calculated self-consistent potential profiles and electron density distributions. Based on these results 
we have shown that - in certain bias voltage and magnetic field ranges - the spin polarization of the 
current exhibits the plateau behavior with the nearly full spin polarization. This property is very 
promising for possible applications in spintronics. 

PACS numbers: 72.25.Dc, 85.30.Mn 
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I. INTRODUCTION 

A rapid progress in homo- and heteroepitaxy of semi- 
conductors with magnetic dopants^"^ has led to the in- 
terest in the transport properties of the diluted magnetic 
semiconductors (DMSs) due to their potential spintronic 
applications. A special attention is paid to the studies 
of (Zn,Mn)Seiaai and (Ga,Mn)AsJ^-i^. The nanostruc- 
tures consisting of the non-magnetic and magnetic lay- 
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lead to a possibility of creation of the nanodevice, 
in which spin polarization of the current can be controlled 
by the bias voltage. The fabrication of an effective spin 
filter for application to semiconductor spintronic nanode- 
vices, such as spin transistor", is one of the most impor- 
tant challenge for the contemporary nanotechnology. 

The application of DMS nanostructure as an effec- 
tive spin filter has been first proposed by Egues^^ft The 
nanostructure studied in Ref.^° consisted of a paramag- 
netic semiconductor layer made of Zni_2:Mn2:Se, sand- 
wiched between two non-magnetic ZnSe layers. In the 
presence of the external magnetic field, the giant Zee- 
man effect occurring in the Zni_a;Mn2:Se layer— causes 
that the paramagnetic layer acts as a potential well for 
spin-down electrons and as a potential barrier for spin- 
up electrons, which leads to the total current dominated 
by the spin-down electrons. In this type of the nanode- 
vice, the spin polarization of the current can be controlled 
by the external magnetic field. There is another type 
of the nanodevice, namely, the resonant tunneling diode 
(RTD) with the paramagnetic quantum well, in which 
the spin polarization of the current can be controlled by 
the bias voltage in the presence of the external magnetic 
field. If the magnetic field is applied to the paramag- 
netic RTD, the resonant tunneling conditions are satis- 



fied for the different bias voltages for the spin-up and 
spin-down electrons due to the giant Zeeman splitting 
of the quasi-bound state in the paramagnetic quantum 
well. This leads to the separation of the corresponding 
resonant current peaks and consequently to the spin po- 
larization of the current. The paramagnetic RTD based 
on the ZnBeSe/ZnMnSe heterostructure was experimen- 
tally demonstrated by Slobodskyy et d\. f^ and theoreti- 
cally investigated by Havu et al.^1 In both the experimen- 
tal and theoretical studies, the separation of the resonant 
current peaks corresponding to the different spin compo- 
nents of the current was found. Nevertheless, in these 
studies, we did not find any hint for a current hysteresis 
and its influence on the spin polarization of the current. 
On the other hand, the nonlinear effects, in particular, 
the current hysteresis (current bistability), can be ob- 
served in the resonant tunneling through the nanodevices 
made from the non-magnetic materials.—"— The first ob- 
servation of the current hysteresis in the nonmagnetic 
(AlGa)As/GaAs/ (AlGa)As heterostructure was reported 
by Goldman and Tsuii2^ The current hysteresis^ occurs 
in the negative differential resistance (NDR) region of 
the current-voltage characteristics and is attributed to 
the accumulation of the electrons in the quantum well 
layer. SoUner^'' proposed another explanation of the 
current bistability based on the oscillations of a circuit 
containing an element that exhibits the NDR. In non- 
magnetic RTDs, the following two types of the current 
hysteresis can be distinguished: (I) the hysteresis occur- 
ring in the bias voltage range below the resonant current 
pea k^^i^^ , (II) the hysteresis occurring in the bias volt- 
age range above the resonant current peakf^i Moreover, 
one can observe both the types of hysteresis on the same 
current- voltage characteristics (double hysteresis) 3^ Re- 
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cently, Dai et al<^ have investigated a possibility of tun- 
ing the current hysteresis in the non-magnetic resonant 
tunnehng diode by changing the geometric parameters of 
the nanodevice and found that the increase of the cohec- 
tor barrier width can enhance the current hysteresis of 
type (I), while the decrease of the collector barrier width 
enhances the type (II) current hysteresis. These results^^ 
indicate that the geometric parameters considerably af- 
fect the nonlinear transport properties of the nanodevice. 

However, in the paramagnetic RTD, the current hys- 
teresis and its influence on the spin polarization of the 
current has not been studied until now. In this paper we 
present the results of such investigations. We show that 
the current hysteresis is of crucial importance for obtain- 
ing the spin-polarized current. We provide the physical 
interpretation of this effect based on the analysis of the 
self-consistent potential profile and the electron density 
distribution. We have predicted the bistable behavior of 
the spin current polarization as a function of the bias 
voltage and the appearance of plateaus of the spin cur- 
rent polarization as a function of the external magnetic 
field. 

The paper is organized as follows: in Sec. 2, we 
describe the model of paramagnetic RTD and the self- 
consistent Wigner-Poisson method. Section 3 contains 
the results, Section 4 - the discussion, and Section 5 - 
the conclusions and summary. 
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FIG. 1. (Color online) Potential energy profile in the para- 
magnetic RTD for the electrons with spin-up (solid line, red) 
and spin-down (dotted line, blue). Coordinate z is measured 
along the layer growth direction, ^l(r) is the electrochemical 
potential of the left (right) contact. The undoped (active) 
region consists of the paramagnetic quantum well made of 
Zni_a;Mna;Se sandwiched between two Zno.95Beo.05Se poten- 
tial barrier layers. The active region of the nanodevice is 
separated from the n-doped ZnSe contacts by the two spacer 
layers at the left and right side. 



II. MODEL AND THEORY 

We consider the paramagnetic RTD diode consisting 
of the Zni_2;Mn2;Se paramagnetic layer sandwiched be- 
tween the two Zno.95Beo.05Se layers (Fig. [T]). The active 
region of the nanodevice is separated from n-doped ZnSe 
contacts by two spacer layers located at the left and right 
contacts. The difference between the conduction band 
minima of Zni_a;Mna;Se and Zno.95Beo.05Se leads to the 
potential energy profile with two barriers and one quan- 
tum well, in which the quasi-bound states can be formed. 
We focus on the effect of bias voltage H applied between 
the left and right contacts and the external magnetic field 
B = (0, 0, B) applied in the growth direction on the spin 
polarization of the current. In the presence of the mag- 
netic field, the exchange interaction between the spins of 
the Mn^+ ions and the spins of the conduction band elec- 
trons leads to the giant Zeeman splitting of the conduc- 
tion band minima in the paramagnetic quantum wellAS. 
This splitting occurs at temperatures below the Curie 
temperature of Zni_2;Mn2,Se and gives rise to the differ- 
ent potential energy profiles for spin-up and spin-down 
electrons. For the small concentration of Mn^+ ions the 
giant Zeeman splitting can be expressed by the formulaic 

where A^oc* = 0.26 eV is the sp-d exchange constant for 
the conduction band electronsjSE x is the concentration 



of Mn^+ ions, SqBs is the effective Brillouin function for 
spin S — 5/2 that corresponds to the spin of Mn^+ ion, g 
is the effective Lande factor, /zb is the Bohr magneton. So 
and Tcff are the phenomenological parameters accounting 
for Mn-Mn interaction that are chosen as foUows:^^ Sq = 
1.60 and T^g = 1-5 K for a; = 0.043, while 5*0 = 1.18 and 
Tcff = 2.55 K for a: = 0.083. 

The potential energy profiles for the spin-up and spin- 
down electrons in the nanodevice are displayed in Fig. [T] 

We calculate the current- voltage characteristics of the 
paramagnetic RTD using the modified version of the 
self-consistent Wigner-Poisson approach. For this pur- 
pose we extend our previous approach^ by assuming the 
two-current modelf^i^ according to which the conduc- 
tion band electrons are described by the spin-dependent 
Wigner distribution function (WDF). In order to cal- 
culate the electronic transport through the nanodevice 
we apply the simplified version of the quantum kinetic 
equation. Assuming the translational invariance in 
the X — y plane, the quantum transport equations for 
the steady state can be reduced to the following one- 
dimensional form: 

m oz znn J 

where p'^{z,k) is the spin-dependent WDF, k is the z- 
component of the wave vector, m is the electron conduc- 
tion band mass, and a =t, 4, is the electron spin index. 
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The non-local Wigner potential for spin channel a is 
given by the formula 



U^{z,k~k')^ / dz' [U^{z + z'/2) 



(3) 



- U„{z - z'/2)] exp [~i{k- k')z'] , 

where Ua-{z) is the total spin-dependent potential energy 
that consists of the two terms: conduction band potential 
energy U^{z,B) and Hartree potential energy U^{z). 

At this stage of our study, we neglect the exchange 
interaction between the conduction band electrons. The 
conduction band potential energy has the form 

N 

U°{z; B) = J2 U^Q{z - z,)e(z,;+i -z)+ (4) 
1=1 

+ A{B)Q{z ~ zi)Q{z2 ~ z) , 

where Zi is the position of the barrier/ well 
(paramagnetic/non-magnetic) interface, 0{z) is the 
Heaviside step function, Ui is the height of the i-th 
barrier and A(i?) is the giant Zeeman splitting of the 
conduction band in the paramagnetic layer. 

In Eq. (jH), the Zeeman splitting for the conduction 
band electrons has been neglected because it is a few or- 
ders of magnitude smaller than the giant Zeeman split- 
ting occurring in the paramagnetic quantum well. 

The Hartree potential energy satisfies the Poisson 
equation 

d^U^iz) 



-[NDiz)~n,{z)] 



(5) 



dz^ EqE 

where Sq is the vacuum electric permittivity, s is the rel- 
ative static electric permittivity, and Nd{z) is the con- 
centration of the ionized donors. 

The density of the electrons with spin a can be ex- 
pressed as follows: 



(6) 



Quantum kinetic equation ([2]) and Poisson equation ([5]) 
constitute a system of equations that can be solved by 
a self-consistent numerical procedure. '^^ We assume the 
Dirichlet boundary conditions for the Poisson equation, 
i.e., f/^(0) = and U^{1) = -eVt, where Vb is the 
bias voltage applied between the left (L) and right (R) 
electrodes separated by distance I. 

When solving quantum kinetic equation ^ we use the 
open boundary conditions'^^ in the form 



P^iO,k) 



(7) 



fc>0 



fc<0 



fRa{k) , 



where f„a{k) {i' = L, R) is the supply function'^^ that for 
contact ly has the form 

f (i.^ t^^bT ^ r r 1 fh'e 

Ju<j[k) = — In S 1 + exp 



ksT \ 2m 



where T is the temperature, and /ijyo- is the electrochem- 
ical potential of reservoir v. 

The current density for non-interacting spin channels 
can be expressed by the following formula: 



e 

2^ 



T)k 

dk—p'^{z,k-Vb,B) . (9) 
m 



(8) 



We note that although the right-hand side of Eq. ^ con- 
tains coordinate z, the current is independent of z for the 
steady state solution. The present calculations have been 
performed with the geometric parameters of the nanode- 
vice, for which the most pronounced current hysteresis 
has been found in non-magnetic RTD 3^ We have checked 
that - in the paramagnetic RTD - the current hysteresis 
effects are also the strongest for the same geometry as 
in nonmagnetic RTD. In the nanodevice studied in the 
present paper, the contacts are made from the n-type 
ZnSe with the homogeneous ionized donor concentration 
Nd = 2 X 10^* cm~'^, the thickness of each contact is 
equal to 17 nm, the thickness of each spacer layer is 3 nm, 
the thickness of each potential barrier (well) layer is 3 nm 
(5 nm), the total length of the nanodevice / = 54 nm, 
and the height of the potential barrier Uq = 0.115 eVj^ 
The energy of the left-contact conduction band bottom 
is taken as the reference energy and set equal to 0. Due 
to the small thickness of the double-barrier region we as- 
sume that the electrons are described by the conduction 
band effective mass of ZnSe, i.e., m — 0.16 toq, where 
mo is the free electron mass. We take on the relative 
electric permittivity e = 8.6 and the lattice constant 
a = 0.5667 nm for ZnSe. The present simulations have 
been carried out for temperature T = 1.2 K. We have ap- 
plied the computational grid with = 95 mesh points 
for coordinate z and = 72 mesh points for wave vector 
k. 



III. RESULTS 

Using the steady state Wigner-Poisson approach, de- 
scribed in Sec. HIl we have performed calculations of 
current-voltage characteristics I{Vh) for the paramag- 
netic RTD in the presence of external magnetic field for 
the spin-up and spin-down current components, sepa- 
rately. Each of these characteristics have been obtained 
for the two directions of the bias sweep: forward bias 
sweep (FBS), in which the bias voltage increased from 
V to 0.2 V with step 0.005 V and backward bias sweep 
(BBS), in which the bias voltage decreased from 0.2 V to 
V with step 0.005 V. In Fig. [51 we present the calcu- 
lated current-voltage characteristics for different values 
of the magnetic fields. 

The results of Fig. [2] show that the resonant current 
peaks corresponding to the spin-up and spin-down cur- 
rent components for FBS and BBS are separated. The 
separation of the resonant current peaks results from the 
giant Zeeman splitting of the energy levels correspond- 
ing to the quasi-bound states formed in the paramagnetic 
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FIG. 2. (Color online) Current-voltage characteristics of the paramagnetic RTD for spin-up (red) and spin-down (blue) currents 
for magnetic field (a) B = 2 T, (b) S = 4 T, (c) B = 6 T, and (d) B = 8 T. The results for the forward (FBS) and backward 
(BBS) bias sweep are displayed by solid line and dotted curves, respectively. The directions of the bias sweep are marked by 
the arrows. In panel (c), the first (I) and the second (II) bias voltage range, in which the current hysteresis occurs, are marked. 



quantum well and increases with the increasing magnetic 
field. This means that the resonance conditions for the 
spin-up and spin-down electron currents are satisfied for 
the different values of the bias voltage. In order to get 
a more deep physical insight into this effect, we present 
in Fig. [3] the spin-dependent electron density for mag- 
netic field B = 6 T in the case of FBS and for the 
following two values of the bias voltage: Vt = 0.075 V 
[Fig. [3] (a)] that corresponds to the current peak for spin- 
down current component and Vb = 0.12 V [Fig. [3] (b)] 
that corresponds to the current peak for spin-up current 
component. The larger accumulation of the spin-down 
electrons in the quantum well [Fig. 3(a)] indicates that 
the resonance condition is better satisfied for these elec- 
trons. In Fig. 3(b), the electron density distributions for 
spin-up and spin-down electrons are completely different 



in the quantum well layer. In this case, we have obtained 
the fairly high concentration of spin-up electrons with 
a simultaneous almost complete depletion of spin-down 
electrons in the quantum well. This electron density dis- 
tribution corresponds to the resonance conditions satis- 
fied for the spin-up electrons only. 

The results plotted in Figs. [2lja-d) show that the in- 
creasing magnetic field shifts the position of the current 
maximum for the spin-up current component towards the 
higher values of the bias voltage, while the position of the 
current maximum for the spin-down current component 
is shifted towards the lower values of the bias voltage. 
Fig. H displays the separation AV = VP^"'' - V^Zn be- 
tween the resonant current peaks for the spin-up and 
spin-down current components as a function of magnetic 
field. The separation between these peaks increases with 
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FIG. 3. (Color online) Electron density distribution for FBS 
at magnetic field B = 6 T. Panel (a) shows the results for 
Vb — 0.075 V that corresponds to the resonant current peak 
for spin-down electrons, panel (b) shows the results for Vt — 
0.12 V that corresponds to the resonant current peak for spin- 
up electrons. 



the increasing magnetic field and reaches ~ 50 mV at 
B = 8 T. For each value of the magnetic field AV for 
BBS is lower than that for FBS. 

The efhciency of the spin filter is characterized by the 
spin polarization of the current defined as follows: 



Jt-Ji 

h + h 



(10) 



where and are the current densities for the spin-up 
and spin-down electrons, respectively. 

Fig.[2]shows that the current- voltage characteristics for 
the spin-up and spin-down current components obtained 
for the increasing bias (FBS) essentially differ from those 
obtained for the decreasing bias (BBS). This leads to the 



FIG. 4. (Color online) Separation AV between the resonant 
current peaks for the spin-up and spin-down current compo- 
nents as a function of magnetic field B in the case of FBS 
(solid, red) and BBS (dashed, blue). 



hysteresis that occurs for each spin component of the cur- 
rent. These hysteresis loops and the phenomena respon- 
sible for them determine the current polarization. The 
calculated spin current polarization is depicted in Fig. [5] 
as a function of bias voltage, and Fig. [5] as a function of 
magnetic field. We observe the distinct plateau regions 
of the spin current polarization that appear for certain 
bias voltage (Fig. [5]) and magnetic field (Fig.lH]). In some 
ranges of the bias voltage and magnetic field, we have ob- 
tained almost full spin polarization of the current flow- 
ing through the paramagnetic RTD. The different spin 
current polarization obtained for FBS and BBS [Fig. [S] 
results from the current hysteresis (cf. Fig. [2]). For ex- 
ample, for Vb = 0.1 V and B = 6 T the spin current 
polarization Pj ~ 1 in the case of the BBS but Pj ~ 0.5 
in the case of the FBS. Fig. [5] shows that the spin cur- 
rent polarization can be large, i.e., 0.8 < Pj < 1, and 
the spin-up current component dominates for both the 
FBS and BBS. In the case of the BBS, the bias voltage 
range corresponding to the almost full spin current po- 
larization is shifted toward the lower voltage. The closer 
inspection of the Pj{Vb) curve for the FBS shows that we 
have obtained the plateau behavior of the spin current 
polarization in quite wide bias voltage ranges [cf. Fig. [5] . 
The width of these ranges increases with the increasing 
magnetic field. In Fig. [5l the bias voltage ranges corre- 
sponding to the plateaus of the spin current polarization 
are marked by the vertical dotted lines. 

The plateau behavior of the spin current polarization 
also occurs if we change the magnetic field keeping the 
bias voltage fixed. This effect is demonstrated on Fig. [51 
which displays the results of calculations for the FBS. 
Depending on the bias voltage we obtain the dominating 
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FIG. 5. (Color online) Spin polarization Pj of the current as 
a function of bias voltage Vb for the FBS and BBS and for 
different values of magnetic field B. The bias voltage ranges 
corresponding to the plateaus of the spin current polarization 
are marked by the vertical dotted lines. 



current of the spin-down electrons for V < Vf, < 0.08 V, 
while for Vb > 0.08 V the spin-up electrons give the main 
contribution to the current. While the spin-down current 
polarization monotonically increases with the increasing 
magnetic field, the spin-up current polarization rapidly 
increases at low magnetic field and becomes almost con- 
stant in the range of the intermediate and high magnetic 
field. In Fig. [6l we can distinguish the following ranges 
of the magnetic field, in which the plateau behavior of 
the spin current polarization occurs: from 1 T to 6 T for 
Vb = 0.1 V and from 0.6 T to 6 T for Vf, = 0.135 V. For 
V{, = 0.115 V we have obtained the three magnetic-field 
ranges with almost constant spin current polarization: 
the first from T to 1.6 T with Pj ~ 0, the second from 
1.6 T to 3 T with Pj ~ 0.5, and the third from 3 T to 
6 T with P, ~ 0.95. 



IV. DISCUSSION 

The calculated current-voltage characteristics [Fig. [2] 
show the hysteresis for each spin component of the cur- 
rent. We have found that the current hysteresis occurs 
in the following two separate bias voltage ranges: (I) the 
first one with the bias voltage values below that the corre- 
sponding to the resonant current peak and (II) the second 
with the bias voltage values above that corresponding to 
the resonant current peak [cf. Fig. (He)]. In these two 
ranges, the current hysteresis originates from two differ- 
ent effects. In region (I), the current hysteresis results 
from the accumulation of the electron charge in the cen- 
tral quantum well, i.e., that made from the Zni-ajMn^jSe 



FIG. 6. (Color online) Spin polarization Pj of the current as 
a function of magnetic field B for the FBS and different bias 
voltage Vb- 



[cf. Fig. [T] . In region (II) , the current hysteresis origi- 
nates from the creation of the quasi-bound state in the 
quantum well induced in the region of the left contact 
and the resonant electron tunneling through this state. 

In order to get a more deep physical insight into the 
origin of the current hysteresis we consider the current- 
voltage characteristics for spin-up component of the cur- 
rent at the magnetic field B = 6T for which the electron 
density distributions and potential energy profiles are cal- 
culated [Fig. [7fa,b)]. First, we discuss the effect of the 
current hysteresis in the first bias voltage range below the 
resonant current peak. We see [Fig. [Tja)] that the poten- 
tial energy profiles and the electron density distributions 
obtained for the increasing and decreasing bias voltage 
are different. In the case of FBS, the electrons are ac- 
cumulated in the paramagnetic quantum well, which is 
in contrast to the case of BBS, for which the electron 
density in the quantum well is low. Now, we consider 
the processes that occur if the bias voltage is gradually 
changed. The increasing bias voltage (FBS case) leads 
to the lowering of the potential well energy, which causes 
that the energy of the quasi-bound state in the param- 
agnetic quantum well also becomes lower. As a result 
at the certain bias voltage, the resonant tunneling con- 
dition is reached. The resonant tunneling appears if the 
energy of the quantum-well quasi-bound state falls into 
the transport window, i.e. it takes the value between the 
minimum of the conduction band for n-ZnSe and the elec- 
trochemical potential of the left contact. In the resonance 
tunneling regime, the current density reaches the maxi- 
mum and we observe the accumulation of the electrons 
in the central quantum- well region [Fig. [7fa)]. The elec- 
trons accumulated in the quantum well create the electric 
field that shifts up the potential well bottom and conse- 
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FIG. 7. (Color online) Electron density distributions and po- 
tential energy profiles for spin-up electrons at magnetic field 
B = 6 T for the bias (a) VJ, = 0.12 V corresponding to range 
(I), (b) Vb = 0.145 V corresponding to range (II) (see text). 
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quently the energy of the quasi-bound state. Therefore, 
the effect of the further increase of the bias voltage that 
leads to getting out the resonance is partly compensated 
by the accumulation of electrons in the quantum well. 
This causes that - in the case of FBS - the resonant tun- 
neling condition is satisfied in a wider bias voltage range. 
As might be expected this effect is not observed for the 
BBS case because the decreasing bias voltage and the 
charge accumulated in the quantum well shifts up the 
potential well bottom, which increases the energy of the 
quasi-bound state in the quantum well region. Therefore, 
the different effects of the increasing and decreasing bias 
voltage on the quantum transport conditions allow us to 
explain the appearance of the current hysteresis in bias 
voltage range (I). 

The current hysteresis in bias voltage range (II) lo- 
cated above the bias corresponding to the resonant cur- 
rent peak manifests itself as the increase of the current 
in the NDR region on the current-voltage characteristics 
[cf. Fig. [5]. In this range of bias voltage, the poten- 
tial energy profiles and electron density distributions ob- 
tained for the increasing and decreasing bias voltage are 



FIG. 8. (a) Expectation value of the energy of the quasi- 
bound state localized in the central quantum well and in the 
quantum well formed in the left contact calculated for the bias 
voltage from region (II) [Fig. 2(c)]. (b) Schematic illustration 
of the resonant tunneling through the quasi-bound states of 
the left contact potential well and central potential well. This 
process is responsible for the current hysteresis in bias voltage 
range (II) above the resonant current peak. 



different (Fig. [71[b)), they also exhibit an essential differ- 
ence in comparison to those obtained for range (I) [cf. 
Fig. El^a)]. In the case of FBS, we observe the depletion 
of the electron density in the region of the left contact 
and the formation of the shallow potential well in this 
region. In the NDR region the system rapidly escapes 
from the resonant tunneling condition and consequently 
the reflection probability from the left barrier drastically 
increases. The wave function of the electron incoming 
from the left reservoir interfere with the wave function of 
the electron reflected from the left barrier, which leads 
to the electron density oscillations in the region of the 
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left contact [cf. Fig. [ZIb)]- In a self-consistent man- 
ner, this electron density distribution leads to creation 
of the shallow potential well, in which the quasi-bound 
state can be formed. The results of calculations of the en- 
ergy expectation values for the quasi-bound states in the 
central quantum well and in the quantum well located 
in the left contact [Fig. [SJa)] allow to explain the rela- 
tion between these quasi-bound states. It appears that 
for each bias voltage value from region (II) the energy 
of the quasi-bound state localized in the region of the 
left contact is lower than the energy of the quasi-bound 
state localized in the central quantum well. If the bias 
voltage increases, the difference between these energies 
decreases from ~ 10 meV for Vt ~ 0.125 V to 2 meV 
for Vb — 0.153 V. We note that this energy difference is 
small and comparable to the energy broadening of the 
corresponding energy level. Therefore, the resonant tun- 
neling conditions are satisfied for some range of the bias 
voltage and the electrons can tunnel via the quasi-bound 
states from the left contact through the central quantum 
well to the right contact leading to the linear increase of 
the current in the NDR region. The resonant tunneling 
that occurs for the FBS in region (II) is schematically 
illustrated in Fig. Hfb). In the case of BBS, the potential 
well in the left contact region is not deep enough in order 
to the quasi-bound state to be formed [cf. Fig. E^b)]. 
Therefore, the tunneling can appear with a very small 
probability and the BBS current is very small. These ef- 
fects lead to the current hysteresis in bias voltage range 
(II) located above the resonant current peak. 

The plateau behavior of the spin current polarization 
that appears in certain ranges of the bias voltage (Fig. 
[51) and magnetic field (Fig. ^ are the most interest- 
ing results of the present calculations. We demonstrate 
that all of these plateaus are caused by the tunneling of 
spin-polarized electrons through the quasi-bound state in 
the region of the left contact. As we mentioned before 
the same process leads to the increase of the current in 
the NDR region and consequently the current hystere- 
sis occurring in bias voltage range (II). For this purpose 
we analyze the spin-dependent electron distributions and 
potential energy profiles calculated for different bias volt- 
age values at the fixed magnetic field (Fig. |9]) and for 
different values of the magnetic fields at the fixed bias 
voltage (Fig. [TU]). Fig. [HI displays the results for the two 
bias voltage values that correspond to the two plateaus 
of spin current polarization in the lower and higher bias 
voltage regions (cf. Fig. [5])- Fig. IHl^a) shows that - for 
the lower bias voltage and for the spin-down electrons - 
we deal with the charge depletion in the central quan- 
tum well and the formation of the quantum well in the 
region of the left contact. In this quantum well, the spin- 
down electrons can be accumulated in the weakly bound 
states. In the higher bias voltage region [Fig. [SJb)], the 
roles of electron spin states are interchanged. As a re- 
sult, the charge depletion appears for the spin-up states 
in the central quantum well and the spin-up electrons 
become accumulated in the left-contact related poten- 
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FIG. 9. (Color online) Electron density distributions and po- 
tential energy profiles for spin-up and spin-down electrons in 
the case of FBS at magnetic field B = 6 T for bias voltage 
(a) Vb = 0.009 V and (b) Vt, = 0.14 V. 



tial well. Therefore, in both the lower and higher bias 
voltage regions, the almost constant spin current polar- 
ization results from the resonant tunneling of the spin- 
polarized electrons via the quasi-bound state formed in 
the region of the left contact. The current of the spin- 
down electrons is responsible for the plateau in the lower 
bias voltage range, while the plateau of the spin current 
polarization in the higher bias voltage range is due to the 
resonant tunneling of the spin-up electrons. 

The spin current polarization versus magnetic field 
dependence exhibits three plateaus for the bias voltage 
Vb = 0.115 V (cf. Fig. [B]). In order to explain the ori- 
gin of this interesting behavior we have calculated the 
spin-dependent electron distributions and potential en- 
ergy profiles for magnetic fields B ~1 T, 2 T, and 4 T 
that correspond to the first, second, and third plateau, re- 
spectively, on the Pj{B) curve (Fig.[TU]). In the magnetic 
field range corresponding to the first plateau fFig.lTUTa)). 
we have obtained the depletion of both the spin-up and 
spin-down electrons in the central quantum-well region. 
Moreover, in the region of the left contact, the shallow 
quantum well is formed for both the spin-up and spin- 
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down electrons. In this left-contact related potential well, 
the quasi-bound states of the electrons with either spin 
can be created. These quasi-bound states participate in 
the resonant tunneling through the paramagnetic RTD 
by means of the same process as observed in case of the 
electron tunneling in the range of the current hysteresis 
marked (II) on Fig. ^c) . 

Since the potential energy profiles and electron den- 
sities for the spin-up and spin-down electrons are nearly 
the same, both the spin currents are almost equal to each 
other, which leads to the vanishing spin current polar- 
ization [cf. Fig. rTOT a)]. The increasing magnetic field 
changes - in different manner - the resonant tunneling 
conditions for each spin polarization of the current. If 
magnetic field exceeds B = 1.6 T, the quantum well for 
spin-up electrons in the region of the left contact disap- 
pears [Fig. [TUTb)]. At high magnetic fields, the spin- up 
electrons become accumulated in the central quantum- 
well layer. This means that in the region of the central 
quantum well the resonant tunneling condition is satis- 
fied only for the spin-up electrons, whereas the spin-down 
electrons can tunnel through the quasi-bound state in 
the region of left contact. Therefore, the spin-up current 
component becomes large, while the spin-down compo- 
nent is only slightly changed (cf. Fig. As a result of 
these changes, we obtain - at the intermediate magnetic 
field (Fig. in]) - the second plateau of the spin current 
polarization with Pj ~ 0.5. If the magnetic field exceeds 

3 T, the left-contact related quantum well also dis- 
appears for the spin-down electrons. At high magnetic 
fields, the density of the spin-up (spin-down) electrons is 
large (small) in the central quantum well. In this case, 
the resonant tunneling condition is satisfied for the spin- 
up electrons, while the spin-down electrons are out of the 
resonance. Therefore, the spin- up current reaches the 
maximum while the spin-down current becomes negligi- 
bly small [cf. Fig. [H^d)]. This leads to the third plateau 
on the Pj [B) curve with ~ 1, i.e., we obtain the almost 
full spin current polarization in the high magnetic field 
ranges (Fig. [S]). We note that for the lower and higher 
bias voltage the system makes the direct transition to the 
fully polarized spin current at certain magnetic field (cf. 
Fig.©. 
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V. CONCLUSIONS AND SUMMARY 

In the present paper, we have studied the spin- 
dependent electron transport through the paramagnetic 
RTD based on Zni_j;Mna;Se/Zno.95Beo.o5Se heterostruc- 
ture. The current-voltage characteristics of the nanode- 
vice have been determined by the numerical solution of 
self-consistent Wigner-Poisson model of electron trans- 
port in the framework of the two-current model. The 
application of the two-current model is fully justified in 
the case of the nanodevice because the magnetic field pro- 
duces the giant Zeeman splitting [Eq. ([T])] which modifies 
the potential energy profile and as a result the electrons 



FIG. 10. (Color online) Electron density distributions and 
potential energy profiles for spin-up and spin-down electrons 
in the case of FBS for bias voltage Vb ~ 0.115 V and magnetic 
fields (a) B = 1 T, (b) B = 2 T, and (c) B = 4 T, which 
correspond to the first, second, and third plateau of the spin 
polarization displayed on Fig. 6. 



with opposite spins independently tunnel via the different 
quantum- well states. Numerical estimation of the ratio of 
the cyclotron resonance energy hujc, where lOc = eB/rUe, 
to the Fermi energy of n-doped ZnSe for maximum value 
of magnetic field {B = 8 T) which is considered in this 
paper gives ^0.46. It means that the relative strength of 
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the maximum magnetic field is rather weak and therefore 
we can neglect the effect of magnetic field quantization. 
Due to the complexity of the numerical calculations, we 
have neglected the intervalley coupling and effects of the 
electron scattering processes which result from the inter- 
action of electrons with phonons and impurities. 

The calculated current-voltage characteristics exhibit 
two separated resonant current peaks that correspond 
to the spin-up and spin-down current components. We 
have found that the current hysteresis appears on the 
current-voltage characteristics for each spin polarization 
of the current. The current hysteresis occurs in two dis- 
tinct bias voltage regions: the first with the bias voltage 
values below that corresponding to the resonant current 
peak and the second - above the resonant current peak. 
The analysis of the self-consistent potential energy pro- 
files and electron density distributions in the nanodevice 
allows us to give the physical interpretation of these ef- 
fects. We have demonstrated that the current hysteresis 
in the first bias voltage region is related to the charge 
accumulation in the central quantum- well layer, while 
the current hysteresis in the second bias voltage region 
results from the resonant tunneling through the quasi- 
bound state created in the region of the left contact. We 
have shown that the current hysteresis leads to the differ- 
ent spin polarization for the forward and backward bias 
sweeps, which - in turn - leads to the appearance of the 
plateaus that appear on the spin current polarization de- 
pendence on the bias voltage and magnetic field. For cer- 
tain plateaus the current is almost fully spin-polarized. 
We have determined the conditions for the full spin cur- 
rent polarization and discussed the underlying physics. 
This property is very promising for future spintronic de- 
vices. 

We expect that including the scattering of electrons 



should result in more smooth current-voltage charac- 
teristics and the decrease of the spin current polariza- 
tion. Moreover, the exchange interaction between the 
conduction-band electrons has been neglected in the 
present calculations. We have also performed the calcula- 
tions with the exchange interaction taken into account by 
the local density approximation. The preliminary results 
show that the incorporation of the interelectron exchange 
changes the present results only slightly, in particular, the 
maximum current polarization becomes slightly smaller 
and the rapid changes on the current-voltage charac- 
teristics become more smooth. We conclude that the 
exchange interaction between the conduction-band elec- 
trons has no significant influence on the present results. 

In summary, we have shown that the electron cur- 
rent flowing through the non-magnetic/magnetic/non- 
magnetic double barrier heterostructure (paramagnetic 
RTD) can exhibit the hysteresis that appears on the 
current-voltage characteristics for the spin-polarized cur- 
rent. Moreover, we have determined the regions of the 
bias voltage and ranges of magnetic field, in which the 
pronounced spin current polarization occurs in the para- 
magnetic RTD. The results of the present paper allow 
us to predict the conditions, under which we obtain the 
spin-polarized current in the double-barrier heterostruc- 
ture with the manganese doped quantum-well layer. 
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